Spatio-temporal evolution and driving factors of new urbanization in central China based on multi-source data

Urbanization is an inevitable outcome of the development of human society to a certain stage, and it is also an irreversible pattern of the concentration degree of human society. Based on multi-source data such as remote sensing images, ecological environment and socio-economic data, the evaluation index system of new urbanization is constructed from multi-dimensions of population, economy, society, space and ecology. To explore the spatio temporal evolution and driving factors of urbanization in 80 prefecture-level cities in central China from 2013 to 2021 by using entropy method, spatial autocorrelation model and geographic detector. The results show that: (1) The level of new urbanization continues to grow, with the average value rising from 0.1562 in 2013 to 0.2557 in 2021, and the regional differences are obvious, forming a circle structure with Wuhan, Zhengzhou and other provincial capitals as the center and weakening radiation to surrounding cities. (2) The agglomeration of ecological urbanization is significant, and the agglomeration trend is gradually enhanced. The high-high agglomeration areas tend to Xinzhou City, most prefecture-level cities in Hubei Province and some prefecture-level cities in Southern Hunan Province, while the low-low agglomeration areas tend to Changzhi City, most prefecture level cities in Henan Province and some prefecture-level cities in Northern Anhui Province. (3) The night light index, total retail sales of consumer goods, investment in fixed assets, proportion of built-up areas and urban economic density are the main driving factors affecting the level of new urbanization. (4) The interaction of driving factors shows double factor enhancement and nonlinear enhancement effects.


Introduction
Urbanization is an important engine of modernization and economic growth, and is regarded as an inevitable trend of human social development [1].China's urbanization has been named as one of the two key factors affecting human development in the 21st century by Stiglitz, winner of Nobel Prize in economics [2].Since the beginning of the 21st century, the urbanization rate of China's permanent residents has increased from 36.22% in 2000 to 64.72% in 2021, which is far lower than the average level of developed and developing countries and has a large space for development [3].In 2014, the state officially issued the National New Urbanization Plan (2014)(2015)(2016)(2017)(2018)(2019)(2020), marking the major transformation of China's urbanization development [4].The 19th National Congress of the Communist Party of China proposed a new urbanization strategy that focuses on promoting people's urbanization and focuses on improving quality [5].The 20th National Congress of the Communist Party of China continued to propose to vigorously promote the construction of new urbanization, which has become an important engine to promote China's economic development and social progress [6].As an important geographical region in China, central China plays a key role in the coordinated development of regional economy and plays an important role in connecting the east with the west [7].In 2016, the National Development and Reform Commission issued the 13th Five-Year Plan for Promoting the Rise of the Central China, which clearly positioned the central China as a key area in the national new urbanization construction.Under the above background, it is of great significance to study the spatio-temporal evolution of new urbanization in central China and its driving factors for optimizing China's urbanization development pattern and promoting regional coordinated and sustainable development.
At present, the definition of urbanization is not completely unified, but there is some consensus on its connotation.Scholars generally believe that urbanization is mainly characterized by population growth, economic development and spatial expansion [8][9][10][11].It is generally believed that urbanization is a regional spatial process, which shows the spatial expansion of urban scale and form, as well as the dynamic evolution of urban functions and functions.In essence, urbanization is a coupling and coordination process between urban activities and supporting environment, and a spatial expression of man-land relationship regional system.In addition to the internal population and economic spatial agglomeration, urbanization also includes material, energy and information exchange with the background environment.Urbanization, as an evolutionary process of economic and social development, is closely related to economic level, ecological environment and social transformation [12].
Studies on urbanization mainly involve urbanization level measurement, development process, influencing factors and urban problems [13][14][15][16].Urbanization level measurement is mainly based on single population data or the construction of evaluation index system to measure urbanization development level [17,18].The measurement criteria of development process are composed of urbanization rate, population density and urban density [19,20].The characteristics and dynamic evolution of urbanization pattern can be revealed through development process [21], and regional differences of urbanization levels at different scales can be analyzed [22].The analysis of influencing factors is mostly carried out from the aspects of economic development [23,24], population migration [25], industrial structure [26], administrative policy and development planning [27].In addition, interdisciplinary integration of urbanization studies with sociology, ecology and other disciplines focuses on the transformation from economic development to sustainable development issues such as resources, environment and rural construction [28,29].This paper explores the sustainable development of new urbanization from the perspectives of the coupling between new urbanization and economic development, ecological environment, basic public services and other factors [30,31], urban diseases caused by urban expansion, ecological resource problems and other development problems and countermeasures [32,33].However, the above studies pay more attention to urban development in terms of population, economy, space and social services, but not enough attention to urban ecological environment, and relatively unified views have not been formed on the selected index system.It often expresses the coordination effect between urbanization and other systems, and cannot accurately express the coordinated development process within urbanization.In addition, most studies focus on the national, provincial and municipal levels, and tend to study coordinated development from a macro perspective, while studies from the middle perspective of prefecture-level cities are rare.
Therefore, based on the panel data of 80 prefecture-level cities in the central China from 2013 to 2021, this paper constructs a comprehensive evaluation index system for new urbanization, uses entropy method to calculate index weights, and uses spatial autocorrelation model and geographic detector to explore the evolution of the spatial-temporal pattern of new urbanization in the central China and its driving factors, so as to provide reference value for the balanced development of the central China.It will provide a scientific basis for the formulation of policies for coordinated development of urbanization in relevant regions.Compared with previous studies, the innovations are as follows: Firstly, based on the mesoscopic perspective, compared with the macro perspective, mesoscopic perspective is more suitable for research on the prefecture-level city level; Secondly, the use of multi-source data, on the basis of traditional social and economic data, the introduction of satellite remote sensing image data to supplement, forming multi-source data to reflect the level of urbanization more three-dimensional; Thirdly, construct a multi-dimensional index system, and build a comprehensive evaluation index system for new-type urbanization from the five dimensions of population, economy, society, space and ecology to make up for the deficiencies of a single dimension.

Study area
The central China is located in the inland hinterland, including 6 neighboring provinces of Hubei, Hunan, Jiangxi, Anhui, Henan and Shanxi as well as 80 prefecture-level cities under their jurisdiction (Fig 1 ), playing an important role in the national regional development pattern [34].The central China has unique geographical location conditions, superior resource endowment, and rich human resources, natural resources, culture and ecological advantages.In recent years, with the deepening of the strategy of the rise of the central China and the development of the Yangtze River Economic Belt, the central China has achieved rapid economic and social development in all aspects, gradually becoming the fourth growth pole of the country's economic growth.By the end of 2021, the central China will have a land area of 1.03 million square kilometers, a permanent population of about 360 million and a GDP of 25 trillion yuan, accounting for 11 percent, 26 percent and 25 percent of the country's total.In the context of accelerated development of social economy and urbanization, it is particularly timely and necessary to take 80 prefecture-level cities in central China as research samples to analyze the spatio-temporal evolution of new urbanization and its driving factors.

Data sources and processing
The data used in this paper mainly include administrative division vector data, ecological environment data, social and economic data, night light index (NLI) remote sensing image data, normalized difference vegetation index (NDVI) remote sensing image data.Among them, the vector data of administrative divisions came from 1:1,000,000 vector data provided by the National Catalogue Service For Geographic Information.The ecological environment and socio-economic data were obtained from the statistical yearbook, ecological environment Bulletin, national economic and social development Statistical Bulletin of provinces and prefecturelevel cities in the central China from 2013 to 2021, and some data were calculated based on the yearbook data.Individual missing data were calculated by interpolation of adjacent years.

Index system construction.
Traditional studies tend to use population or economic attributes to measure the level of urbanization.From the perspective of the development of cities and towns, the development of cities and towns cannot be separated from the functions of population, space, economy and social services.While new urbanization should pay attention to the development of cities and towns themselves, it should also pay attention to the coordinated development of cities and ecology.Considering the availability, completeness and scientificity of index data, referring to relevant research results [35,36], NLI and NDVI data were introduced on the basis of previous studies, and 25 indicators were selected from the five dimensions of population, economy, society, space and ecology in combination with regional differences and the actual situation of the central China to build a new urbanization level evaluation index system (Table 1).Among them, population urbanization emphasizes the differences in population size and urbanization level among prefecture-level cities, economic urbanization highlights the differences in economic development, industrialization and per capita level among prefecture-level cities, social urbanization focuses on local fixed asset investment, social consumption and budget income levels, and spatial urbanization represents the differences in land development degree among prefecture-level cities. Ecological urbanization focuses on the difference of natural capital and ecological environment quality among prefecture-level cities.In order to eliminate the impact of population size differences among prefecture-level cities, each index adopts per capita index.

Entropy method.
The entropy method is a relatively objective multi-index weighting method.Compared with the subjective weighting method, it can better avoid the deviation caused by human factors and more accurately reflect the importance of indicators [37].Since the meaning and dimension of each index are different, in order to make the index data comparable, the range method is used for dimensionless processing, and then the entropy method is used to calculate the weight.The calculation formula is as follows: Negative indicators : Where, Z ij represents the dimensionless value of the j index X ij of the i prefecture level city, and X min and X max are the minimum and maximum value of the j index respectively.
The proportion of the j index of the i prefecture-level city P ij : Entropy value e j of the j index: Where, e j is the entropy value of the j index, 0�e j �1; k is a constant; m is the number of prefecture-level cities and n is the number of indicators.
Then, the information entropy redundancy g j of the j index: Weight w j of the j index : 2.3.3Spatial autocorrelation model.Spatial autocorrelation is the expression of spatial dependence and heterogeneity generated by geographical location or adjacency, and it is a statistical measure of the spatial distribution structure of various elements in the regional system [38,39].Spatial autocorrelation models are divided into global spatial autocorrelation and local spatial autocorrelation.Global spatial autocorrelation is used to test whether research objects are clustered in space, and local spatial autocorrelation is used to detect spatial clustering degree.
Global spatial autocorrelation is applied to measure whether the new urbanization level in central China has significant spatial correlation, which is usually expressed by Global Moran's I, and the calculation formula is as follows: Where, n is the number of prefecture-level cities; W ij is the spatial weight matrix, which is 1 if adjacent and 0 if not adjacent.The x i and x j are the attribute values of new urbanization in central China on corresponding spatial units respectively, � x i and � x j are the attribute average value, and S 2 is the variance of new urbanization.The value range of Global Moran's I is [1,1].If the value is positive, it indicates that the new urbanization in central China presents spatial clustering.If the value is negative, it indicates a discrete distribution.If the value is 0, it is a random distribution.
Local spatial autocorrelation is adopted to identify the spatial dependence and heterogeneity of new urbanization in central China [40], which is usually represented by Local Moran I.The calculation formula is as follows: If the value I is significantly positive, it indicates that the new urbanization level in central China presents a high-value agglomeration area.If it is significantly negative, the low value cluster area will appear.In addition, the Local Indicators of Spatial Association (LISA) cluster map is a method used for spatial data analysis, whose purpose is to identify spatial clustering patterns.According to the cluster maps of Local Moran I and LISA, the new urbanization in central China can be divided into four spatial correlation types: high-high (H-H) cluster, highlow (H-L) cluster, low-high (L-H) cluster and low-low (L-L) cluster [41].
2.3.4Geographic detector model.Geographical detector is a statistical method to detect spatial heterogeneity and reveal the driving force behind it.It is composed of risk detector, factor detector, ecological detector and interactive detector.Compared with traditional statistical methods, geographical detector can overcome the influence of multicollinearity of independent variables and solve the spatial dependence and heterogeneity caused by scale change [42,43].The calculation formula is as follows: Where, q represents the explanation degree of driving factors to the level of new urbanization, q2[0,1], the greater the value, the greater the impact of the factor on urbanization, and vice versa.L represents the number of layers of the driving factor X; N and σ 2 respectively represent the number of units and variance of the research object, N i and σ i 2 respectively represent the number of spatial units and variance in layer i.
According to the action type, interactive detection can be divided into five categories:nonlinear weakening, single factor nonlinear weakening, double factor enhancement, independent and nonlinear enhancement [44].The interaction types are shown in Table 2.As can be seen from Fig 2, the overall urbanization level in central China continues to grow.The lowest value of the comprehensive index is 12.4969 and the highest value is 20.4565.Compared with 2013, the urbanization level in 2021 increases by 63.7%.From the perspective of all dimensions, except the social urbanization level fluctuates greatly, showing an "updown-up" trend, the urbanization level of other dimensions shows a certain degree of stable growth.Among them, the economic urbanization level increased the most, from 2.0403 in 2013 to 4.8659 in 2021, the level of ecological urbanization is relatively stable.

Analysis of spatial distribution characteristics.
According to the results of the comprehensive index of new urbanization level in central China, the natural break point method of ArcGIS10.7 was applied to divide it into five grades: high, higher, medium, lower and low (Fig 3) [45].In 2013, the level of new urbanization was mainly low, with an average value of 0.1562.The number of prefecture-level cities at high level only accounted for 2.5%, with only Wuhan and Zhengzhou.The number of prefecture level cities at a higher level accounts for 5%, mainly distributed in provincial capitals such as Hefei, Taiyuan, Changsha and Nanchang (Fig 3A).In 2017, the overall level of new urbanization was improved to a certain extent, mainly low level and lower level, with an average value of 0.2034.The number of medium level and high level prefecture-level cities increased (Fig 3B).In 2021, the overall level of new urbanization will be further improved, with the average value of 0.2557 and the number of high level prefecture-level cities increasing to 7.However, Chizhou, Yuncheng, Jinzhong, Xinzhou, Linfen, Lvliang, Zhangjiajie and Suizhou have always been at a low level (Fig 3C).In general, the level of new urbanization in central China showed obvious regional differences during the study period, showing a circle structure with the capital cities such as Wuhan and Zhengzhou as the center and weakening radiation to surrounding cities.

Global spatial autocorrelation analysis.
In order to deeply explore the evolution characteristics of the spatio-temporal pattern of new urbanization in central China, the Arc-GIS10.7 software is used to analyze the spatial correlation and difference of regional urbanization through the Global spatial autocorrelation method, and the Global Moran's I and test results of new urbanization in central China are calculated (Table 3).As can be seen from Table 3, the P values of population, economy, society, space and overall urbanization are almost all greater than 0.05, failing the 95% significance test, indicating that the spatial structure of new urbanization in central China tends to be random distribution in these dimensions, and the regional division of labor and cooperation are not obvious, and the spatial agglomeration is not significant.The Moreland index of ecological urbanization is positive and shows an upward trend through 95% significance test, indicating that there is a significant spatial agglomeration of ecological urbanization, that is, prefecture level cities with higher ecological urbanization tend to cluster relatively, while prefecture level cities with lower ecological urbanization tend to cluster relatively, and the agglomeration trend is gradually increasing.This is due to the persistence of ecological civilization construction in central China in the past decade.Therefore, it is necessary to further analyze the spatial dependence and heterogeneity characteristics of ecological urbanization.From 2013 to 2021, the spatial pattern of ecological urbanization in central China showed obvious agglomeration characteristics of high value and low value.The high-high agglomeration areas are expanding from Wuhan, Xinzhou and Lvliang to Xinzhou, most prefecture-level cities in Hubei and some prefecture level cities in southern Hunan.The spatial difference of ecological urbanization level is small, and the ecological urbanization level of itself and neighboring prefecture-level cities is higher, showing a significant positive correlation.The low-high agglomeration areas are scattered and have large spatial differences.Due to their low ecological environment quality, the development of ecological urbanization is limited, resulting in lower ecological urbanization than neighboring prefecture-level cities.The high-low concentration areas are scattered, and the ecological urbanization level is higher than that of the neighboring prefecture-level cities due to the rational development of ecological resources, and there is a significant negative correlation with the surrounding prefecture-level cities. Low-low agglomeration areas continue to expand, and eventually tend to Changzhi, most prefecture-level cities in Henan and some prefecture-level cities in northern Anhui.Due to the combined effect of backward development and adverse natural conditions, regional ecological urbanization lacks power, and the ecological urbanization level of its own and neighboring prefecture-level cities is low.

Analysis of driving factors of new urbanization
By analyzing the spatio-temporal evolution characteristics of new urbanization in central China, it is found that the spatial dependence and heterogeneity of new urbanization are not significant, and it is difficult to accurately identify its driving factors by traditional statistical methods.As an effective tool to identify spatial heterogeneity features, geographic detectors can better analyze the causality generated by spatial correlation, diagnose and solve spatial heterogeneity and its causes [46].This paper mainly uses factor detector and interaction detector to explore the influence degree of each driving factor on new urbanization in central China and the influence degree of each factor interaction.

Factor detection analysis.
The driving factors of the new urbanization system are classified and sorted by k-means clustering method [47].With the urbanization level as the dependent variable, population, economic, social, spatial and ecological urbanization as the independent variables, combined with the geographic detector model, the driving factors of the new urbanization level in central China are detected and analyzed by Geodetector software.The results are shown in Table 4.
As can be seen from Table 4, during the study period, population, economic, social and spatial urbanization almost all passed the 95% significance test, while only ecological urbanization failed the test.From the perspective of the development of cities and towns themselves, the development of urbanization cannot be separated from the functions of population, space, economy and social services, although the new urbanization also focuses on the coordinated development of cities and ecology.However, China's ecological civilization construction  started late, and the construction of ecological urbanization in central China has only been more than ten years, so the driving effect of ecological urbanization on the new urbanization is weak.The q value of urbanization in each dimension has gradually changed from population urbanization > spatial urbanization > economic urbanization > social urbanization in 2013 to economic urbanization > social urbanization > spatial urbanization > population urbanization in 2021.This indicates that the focus of determining the level of new urbanization has shifted from the basic level of population and spatial urbanization to the high-level and advanced level of economic and social urbanization.Further, the specific driving factors and their interaction relationship were explored for 25 secondary indexes, and the results were shown in Table 5.In 2013, the q value of X 15 ranked first among all index factors, indicating that the proportion of built-up areas was the main driving force for promoting urbanization development.The q values of X 16 , X 9 , X 11 and X 12 exceed 0.7 and pass the 95% significance test, indicating that urban economic density, night light index, total retail sales of consumer goods and investment in fixed assets have important driving effects on urbanization.The value of other index factors q was less than 0.7, and part of them failed the significance test, showing weak explanatory ability.In 2017, the q value of urban economic density X 16 jumped to the first place.The q values of X 1 , X 5 , X 9 , X 10 , X 11 , X 12 and X 15 are greater than 0.7, and pass the 95% significance test, indicating that the enhancement of economic and social urbanization has promoted the improvement of urbanization level.The explanatory ability of other index factors is weak.In 2021, the q value of X 16 of urban economic density still ranks first.The q values of X 1 , X 9 , X 10 , X 11 , X 12 and X 15 are greater than 0.7, and pass the 95% significance test, indicating that economic and social urbanization is an important driving force to promote urbanization development.The explanatory ability of other index factors is weak.In general, during the study period, night light index X 9 , total retail sales of consumer goods X 11 , investment in fixed assets X 12 , proportion of built-up areas X 15 and urban economic density X 16 are the main driving factors affecting new urbanization in central China.

Interaction detection analysis.
In order to further explore the interaction relationship among driving factors of new urbanization in central China, the index factors that failed the test were removed at the significance level of 5%, and the interaction detector was used to analyze the influence degree of the interaction of driving factors on new urbanization (Fig 6).In 2013, the interaction between X 12 and X 15 was the strongest, with q value of 0.954.X 11 and X 15 , X 10 and X 15 , X 1 and X 15 , and X 9 and X 24 had strong interaction, and q values were all above 0.9, indicating that the influence of dual-factor interaction on urbanization level was enhanced.In 2017, the interaction between X 12 and X 15 was the strongest, with q value of 0.959.X 5 and X 11 , X 11 and X 15 , X 3 and X 11 are stronger, indicating that their synergistic effect is conducive to the promotion of urbanization level.In 2021, X 12 and X 15 have the strongest influence with q value of 0.957, while X 9 and X 12 , X 11 and X 15 , X 1 and X 15 , X 10 and X 17 , and X 5 and X 10 are stronger, indicating that their synergistic effect is an important driving force for urbanization development.To sum up, the interaction of driving factors during the study period presents double factor enhancement and nonlinear enhancement effects, and the q values of the interaction of any two factors on new urbanization are significantly greater than the single factor independent effect, indicating that the double factor interaction will enhance the interpretation of urbanization level.a.2013 b.2017 c.2021

Conclusion and discussion
Based on the mesoscopic perspective of prefecture-level cities and multi-source data such as remote sensing images, ecological environment and social economy, this paper constructs a new urbanization evaluation index system in multiple dimensions, uses entropy method to calculate index weights, and adopts spatial autocorrelation model and geographic detector to explore the spatio-temporal evolution characteristics and driving factors of urbanization level in 80 prefecture-level cities in central China from 2013 to 2021.The main conclusions are as follows: 1.There are obvious regional differences in urbanization, showing a circular structure with Wuhan, Zhengzhou and other provincial capitals as the center and a weakening radiation to surrounding cities.During the study period, the new urbanization level in central China continued to rise, with a growth rate of 63.7%.Social urbanization showed an "upward downward-upward" trend, and economic urbanization had the largest growth rate.The average level of regional urbanization increased from 0.1562 in 2013 to 0.2557 in 2021, generally at a medium and lower level, while Chizhou, Yuncheng, Jinzhong, Xinzhou, Linfen, Lvliang, Zhangjiajie and Suizhou were always at a low level.
2. The agglomeration of ecological urbanization is significant, and the agglomeration trend is gradually enhanced, the spatial pattern of high value and low value agglomeration is obvious.During the study period, the Moran's I index of ecological urbanization in central China increased from 0.2119 in 2013 to 0.3539 in 2021, while the spatial structure of urbanization in other dimensions tended to be random distribution, and agglomeration was not significant.The high-low and low-high ecological urbanization agglomeration areas are scattered.High-high agglomeration areas are expanding, mainly distributed in Xinzhou, most prefecture-level cities in Hubei Province and some prefecture-level cities in southern Hunan Province.Low-low agglomeration areas tend to Changzhi, most prefecture-level cities in Henan province and some prefecture-level cities in northern Anhui Province.
3. Night light index, total retail sales of consumer goods, investment in fixed assets, proportion of built-up areas and urban economic density are the main driving factors affecting urbanization level.During the study period, the q values of night light index, total retail sales of consumer goods, investment in fixed assets, proportion of built-up areas and urban economic density are all greater than 0.7.The focus of determining the level of new urbanization in central China has shifted from the population and spatial levels to the economic and social levels.The enhancement of economic and social urbanization has promoted the improvement of new urbanization level.
4. The interaction of driving factors showed double factor enhancement and nonlinear enhancement effects.During the study period, the interaction of two factors will enhance the interpretation of the level of new urbanization, and the q value of any interaction of two factors is greater than the independent effect of single factor, and the synergistic effect will help promote the improvement of the level of new urbanization.
The research results of this paper are helpful to clarify the evolution trend and driving mechanism of urbanization in central China, and provide reference value for the balanced development of the central China.At the same time, the article still has some shortcomings.
1.In terms of research scale, the analysis is mainly based on the mesoscopic perspective of prefecture-level cities, and there is no more detailed research on the level of urbanization at the micro-scale county level, which needs to be improved in the future.
2. In terms of index selection, Multi-source data such as remote sensing images, ecological environment and social economy data are combined for selection, but the uniformity and representativity of the indicators need to be further enriched and improved.
3. In terms of research content, different urbanization levels of prefecture-level cities can be further classified and zoned in the future, so as to better provide systematic and targeted reference information for the new urbanization development of different prefecture-level cities, and thus improve the application value of the research results.
NLI and NDVI data are derived from Resource and Environment Science and Data Center of Chinese Academy of Sciences.Based on the NPP/VIIRS satellite noctiluminosity remote sensing image data from 2013 to 2021, the annual degree NLI data set of central China is generated with a spatial resolution of 0.004 degrees.Obviously, the intensity of light is positively correlated with socio-economic factors.This quantitative relationship can be used to study the estimation of socio-economic factors.Based on SPOT/VEGETATION satellite remote sensing image data, the annual NDVI data set of central China is generated by the maximum synthesis method on the basis of monthly data with a spatial resolution of 1km.The data set effectively reflected the distribution and change of vegetation cover in central China on both spatial and temporal scales.It has very important reference significance for the research of ecological environment related fields.

Fig 1 .
Fig 1. Location overview and remote sensing image of central China.Note: Based on the standard map No. GS (2016) 2556 downloaded from the National Catalogue Service For Geographic Information, the base map has not been modified.https://doi.org/10.1371/journal.pone.0298099.g001

3. 1 . 1
Analysis of temporal evolution characteristics.According to the evaluation index system of new urbanization level in central China (Table1),the weight is multiplied by the index standardization results to obtain the comprehensive index of urbanization level in central China from 2013 to 2021, as shown in Fig 2.

Fig 3 .
Fig 3. Spatial distribution of new urbanization in central China.Note: Based on the standard map No. GS (2016) 2556 downloaded from the National Catalogue Service For Geographic Information, the base map has not been modified.https://doi.org/10.1371/journal.pone.0298099.g003

3 . 2 . 2
Local spatial autocorrelation analysis.GeoDa095 software was used to calculate the Local Moran's I value of ecological urbanization in central China from 2013 to 2021, and Moran's I scatter diagram was drawn (Fig 4).From the perspective of local spatial autocorrelation, Moran's I values of ecological urbanization were all positive during the study period, and scattered data were concentrated in the first and third quadrants, indicating that the new urbanization in central China formed a relatively stable spatial pattern in terms of ecology.a.2013 b.2017 c.2021.In order to more accurately measure the temporal and spatial evolution characteristics of ecological urbanization, the LISA agglomeration of ecological urbanization was drawn (Fig 5).

Fig 5 .
Fig 5. LISA cluster of ecological urbanization in central China.Note: Based on the standard map No. GS (2016) 2556 downloaded from the National Catalogue Service For Geographic Information, the base map has not been modified.https://doi.org/10.1371/journal.pone.0298099.g005